Dynamic Light Scattering from Semidilute Actin Solutions: 
A Study of Hydrodynamic Screening, Filament Bending Stiffness 
and the Effect of Tropomyosin/Troponin-Binding 
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Quasi-elastic light scattering (QELS) is applied to investigate the effect of the 
tropomyosin/troponin complex (Tm/Tn) on the stiffness of actin filaments. The importance of 
hydrodynamic screening in semidilute solutions is demonstrated. A new concentration dependent 
expression for the dynamic structure factor g{k,t) of semifiexible polymers in semidilute solutions 
is used to analyze the experimental QELS data. A concentration independent value for the bending 
modulus K is thus obtained. It increases by 50% as a consequence of Tm/Tn binding in a 7:1:1 molar 
ratio of actin/Tm/Tn. In addition a new expression for the initial slope of the dynamic structure 
factor of a semifiexible polymer is used to determine the effective hydrodynamic diameter of the 
actin filament. Our results confirm the general relevance of the concept of (intrinsic) semifiexibility 
to polymer dynamics. 
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I. INTRODUCTION 



The main motivation of the present work was to study 
the effect of tropomyosin/troponin (Tm/Tn) binding on 
the bending stiffness of actin filaments. The semifiexible 
filaments with a length of 400 (??) bind along the groove 
between the two twisted strands of the actin filament 
(Figure |^). The binding is influenced by the presence of 
Ca+"'". This effect is of great biological significance, since 
it plays a central role in the regulation of the myosin- 
actin coupling in muscles. 
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FIG. 1. Model of Ca++ mediated Tm/Tn binding to actin 
filaments. Tm/Tn is known to bind along the groove of the 
actin filaments. 



Actin, when polymerized in vitro, forms semifiexible 
macromolecules of contour lengths L up to some 30 /im 
and with a persistence length Lp of at least some /im. 
The system we are interested in is a semidilute solution 
of these macromolecules, i.e. an entangled network, where 
the single actin filaments interact strongly but the free 
volume is still much larger than the excluded volume. 
As a consequence of the great extension and large per- 
sistence length of the molecules, as compared to their 
lateral diameter a (some nm), the semidilute regime is 
unusually large. We probe the samples with scattering 



wavelengths A = 2n/k that are somewhat shorter than 
the mesh size ^ of the network, which, in turn, is shorter 
than the persistence length Lp; i.e. the condition 



(1.1) 



is fulfilled for all of our samples. An electron micrograph 
of a typical system under study is reproduced in Figure 0. 




FIG. 2. Electron micrograph of a 0.4^ actin solution poly- 
merized in vitro. The bar indicates the length of 1 nm. We 
probe the samples with scattering wavelengths \ = 2Ti/k that 
are somewhat shorter than the mesh size ^ of the network, 
which, in turn, is shorter than the persistence length Lp. The 
global motion of all polymers is strongly hindered by the en- 
tanglement. 

The bending stiffness of a filament-like macromolecule 
may be inferred from an analysis of its conformational dy- 
namics. Therefore various experimental techniques have 
been applied to investigate the dynamical properties of 



1 



actin [Q. Solutions and gels were probed by high sen- 
sitivity rheology using torsional |^ and magnetic bead 
rheometers Kas et al. Q and others [g| analyzed 
single labeled filaments by microfluorescence microscopy 
combined with dynamic image processing. There have 
also been several attempts to establish QELS as a quan- 
titative method to probe the static and dynamical prop- 
erties of semiflexible macromolecules (for critical discus- 
sions of the literature see e.g. Refs. |q-01). But in con- 
trast to the dynamic properties of flexible polymers the 
dynamics of semiflexible polymers is not yet very well un- 
derstood. This is mainly a consequence of the difficulties 
caused by the rigid constraint of a (virtually) unextensi- 
ble contour. Models that try to represent the unextensi- 
ble contour honestly have to deal with considerable 
technical difficulties. On the other hand, models that 
relax the constraint like the various modified versions of 
the so called Harris-Hearst-Beals-model ||l^ (for a short 
summary see [|l^) allow for artificial stretching modes 
and predict a gaussian probability distribution for the 
spacial distances of the contour elements. As the most 
promising model with relaxed constraint we regard the 
model described in Refs. | |ll| , p^ ,p[ . We will neither re- 
lax the constraint, nor do we need all the machinery 
developed by Aragon and Pecora in Refs. (How- 
ever, one of our results, Eq. (2_^), may be obtained by 
their method [^.) In the special case, we are interested 
in, wc profit from some simplificatio ns br ought about by 
the separation of length scales, Eq. (1.1), which in turn 
gives rise to a hierarchy of time scales: QELS measures 
the temporal decay of configurational correlations of the 
filaments. For long semiflexible actin filaments in en- 
tangled networks the internal configurational dynamics 
dominates over the center of mass and rotational mo- 
tion of the molecules in time intervals typically probed 
by QELS. In addition, we can restrict ourselves to the 
weakly bending rod limit, where one can approximately 
take the undulations to be transverse to the mean con- 
tour [p"3|J^], which is virtually fixed for the relevant time 
scales. The decay rate of the dynamic structure factor 
becomes thus accessible to a simple physical interpreta- 
tion in terms of the bending modes only, i.e. the local 
bending modes may be studied "in isolation" . This is a 
great practical advantage compared to the general case. 

But the time decay of the structure factor also de- 
pends on the strength of hydrodynamic correlations. To 
obtain quantitative results it is therefore necessary to 
take into account the screening of the hydrodynamic self- 
interaction of the filaments in semidilute solutions (or 
their fi nite l ength in the dilute case). As explained in 
Section II A this is most simply achieved by means of 
a concentration dependent renormalized friction coeffi- 
cient, C_L- 

In contrast to the expression for the dynamic struc- 
ture factor for intermediate times, which takes a simple 
form only after some approximations, its initial slope, re- 
sulting from the quasi-free Brownian fluctuations about 
the equilibrium configuration may be calculated exactly 



if Eq. (1.1) holds. The result is co mpared with previous 
experimental data in Section II B. The good agreement 



of the theoretical prediction with the available data con- 
firms our hypothesis that the intrinsic rodlike structure 
of polymers is a significant detail, which influences the 
dynamic properties. Moreover, the microscopic hydro- 
dynamic diameter a of the filament, which enters the 
dynamic structure factor as a parameter, also appears in 
the expression for the initial decay rate and may thus be 
determined by a measurement of the initial slope. Sec- 
tion |ll^ lists materials and methods, and the new exper- 
imental results are presented in Section |v| . Finally we 
discuss possible refinements to our approach in Section 



II. THEORETICAL BACKGROUND OF QELS 
FROM SEMIFLEXIBLE POLYMERS 

A. Stretched Exponential Decay of the Dynamic 
Structure Factor 

The dynamic structure factor of a chain of length 
L — NAs with N segments located at r„ (n = 1, . . . , N) 
is defined by 



5^^' ^) = ^ ^(exp[zfc(r„(t) - r„,(0))]) . 



(2.1) 



The brackets (. . .) denote the ensemble average over all 
chain conformations, and k is the scattering wave vector. 
On the typical length scales probed by QELS (0.1 ^ 1 
/ini) the conformational dynamics of actin filaments are 
dominated by the bending undulations. Consequently, 
we take the bending energy, given by the contour integral 
Jj^ ds over the local curvature multiplied by the bending 
modulus K, 



(2.2) 



to be the only relevant energetic term in a canonical de- 
scription. This implies that the static mean square end- 
to-end-distance is given by the well known Kratky-Porod 
formula (lij with the persistence length Lp — n/kBT. 
The time decay of the structure factor g{k,t)/ g{k,0) is 
not easily calculated in the most general case. However, 
the particular case under study allows for simplifications, 
which enable us to adapt some of the ideas common in 
the theory of flexible polymers ||l^,|l^ . We start from the 
Langevin equation for a single polymer 



d_ 

dt 



rs{t) 



(2.3) 



where is the stochastic force (white noise) and H± is 
an effective mobility matrix, which takes into account the 
solvent-mediated self- interaction of the filament. This 
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can be understood in analogy to the usual Oseen tensor. 
However, the effective reduction of the degrees of free- 
dom due to the rigid constraint of a fixed contour length 
requires that the local longitudinal motion be projected 
out. In the weakly bending rod limit this is practically 
achieved by a suitable choice of coordinates: longitudinal 
distances are kept fixed, and the bending undulations are 
described by transverse coordinates (cf. Eq. ( |2.4|) below). 
This is how one would address the problem in analogy to 
the classical Zimm model for flexible polymers, and how 
it was attempted previously in Ref. . We shall neglect 
the center of mass and rotational motion, which are slow 
compared to the internal dynamics of t he m olecule as a 
consequence of the scale separation Eq. (lA), the entan- 
glement and hydrodynamic screening. 





Flexible Filaments 



Semiflexible 
Filaments 



FIG. 3. Flexible filaments interact mostly with themselves, 
whereas semiflexible filaments interact with each other even 
at quite low concentrations. The perturbation of a single 
semiflexible filament by its surrounding may be modeled by 
a screening of the hydrodynamic self-interaction. 

For the following it is important to realize that there 
is a profound difference between the hydrodynamic in- 
teraction of flexible polymers and semiflexible polymers 
in semidilute solutions. A flexible polymer is coiled and 
is thus far more likely to interact with itself than with 
surrounding polymers. On the other hand, semiflexible 
filaments like actin are much more stretched and strongly 
interact with each other even at quite low concentrations 
(Figure As the scattering vector k in our experiments 
is large enough in magnitude to resolve single actin fila- 
ments (see Eq. (I.l)), we are interested in the dynamics 
of a single filament and use a mean field approximation to 
model the hydrodynamic interaction with the surround- 
ing, i.e. we introduce a screening of the hydrodynamic 
self-interaction along a single polymer. This is achieved 
by use of the (preaveraged) screened transverse mobility 
matrix. 



-r/A 



STTrjr 



\r){r\ 



(2.4) 



where r := — r^', A is the screening length and rj 
the solvent viscosity. Beyond this length A, the hydro- 



dynamic self-interaction of a filament is weak and corre- 
lations decay rapidly. The explicit form of the projec- 
tor is a consequence of modelling the actin filament as a 
straight rod on length scales smaller than A with respect 
to the hydrodynamic interaction. This should be a good 
approximation for ^ Lp. 

To simplify the calculation we use a kinetic coefficient 



ia/L)- 



log A/a 
47rr/ 



(2.5) 



in place of the tensor H±. (± is a renormalized fric- 
tion, obtained by taking the terms in parentheses in Eq. 
( ^.3[ ) out of the integral and averaging over all segment 
positions s for a rigid rod of diameter a. Replacing the 
mobility matrix by a simple coefficient amounts to setting 
the effective friction for all modes equal to the friction of 
the dominant mode of wavelength A. In this approxima- 
tion interactions between different modes are neglected 
and the mobility for the very short wavelength modes is 
supposedly slightly overestimated. The latter should not 
profit as much by correlations over distances larger than 
their wavelength, as is implied by Eq. (2.4). This ap- 
proximation will allow for a simple expression [Eq. (2.S) 
below] for the structure factor and at the same time cap- 
tures the main effect of the hydrodynamic self-interaction 
of the single polymers as well as their mutual interac- 
tion. A theoretical basis for the above ansatz was already 
worked out by Muthukumar and Edwards and will 
be explained in more detail in a forthcoming paper [po[ . 

A, a and the persistence length Lp are the three char- 
acteristic length scales in terms of which the dynamics 
of a semiflexible polymer in semidilute solution is char- 
acterized. The microscopic cutoff parameter a takes care 
of the finite thickness of the filament. Its value may be 
determined experimentally by the QELS method as de- 
scribed in the next subsection. We suppose that the hy- 
drodynamic screening length A may be identified (up to 
a numerical factor) with the mesh size ^ of the actin net- 
work. 



A^^ 



(2.6) 



(Basically the same relation A ^ c^^l"^ can be derived 
within the effec tive medium approach for rods jl^.) It 
is through Eq. (2^) that the actin concentration Ca ul- 
timately enters Eq. (E]q) for the structure factor. The 



scaling law for the mesh size ^ in Eq. ( p. 6]) should be 
valid in the semidilute regime, when ^ <Lp and the solu- 
tion appears as a random network of almost rodlikc seg- 
ments. This was indeed confirmed experimentally p9|. 
The dimensionless quantity A/a, which could be called 
the "hydrodynamic aspect ratio" of the polymer, is ac- 
tually the only characteristic parameter of the filament 
entering expression ( |2.8| ) for the structure factor besides 
the persistence length Lp (cf. Eq. (^) below). 
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Detailed calculations [E0[ show that the relaxation time 

of the p*'' mode is Tp — \wp) ' '^'^^ would guess 
from dimensional analysis, and that in the intermediate 
time regime^, 



(fcLp)-4/3^^i ^^^^^^ 



(2.7) 



the structure factor factor is to a good approximation 
given by a stretched exponential 



g{k, t) = g{k, 0) exp ( -fflM(^,t)3/4 



The decay rate, 



Ik 



9tt 



,„8/3 r~l/3 



(2.8) 



(2.9) 



depends on the filament stiffness through the persistence 
length Lp and on the actin concentration Ca (or the mesh 
size) through the renormalized friction coefficient It 
should be comp ared with the initial decay rate given be- 
low in Eq. ( 2.11 ) for the dilute case. Both are expected to 
play an important role in many dynamical problems, i.e. 
they are supposed to constitute the characteristic time 
scales of motion. 

Some remarks have to be made on the derivation and 
the domain of validity of Eq. ( p.8[ ) . First of all it is im- 
portant to realize that the local motion of the monomers 
is neither strictly isotropic (like in the random coil limit) 
nor strictly tra nsve rse to the mean contour (as was as- 
sumed in Eq. (^.4|)). The averaging over the different 
orientations can be performed in both limits |^^. Upon 
neglecting the anisotropy of the local segment motion, 
one obtains the simple stretched expo nential from of the 
dynamic structure factor, Eq. (2^). In the opposite 
limit - i.e. for local segment motion strictly transverse 
with respect to the mean contour - averaging over the 
different orientations is performed after calculating the 
anisotropic dynamic structure factor for an individual 
macromolecule, whose mean orientation and center of 
mass position is assumed to be fixed in space and time. 
In our case this is realized by the time scale separation 
mentioned above and by the fact that the macromolecules 
are imbedded in a network. It is found that the resulting 
dynamic structure factor is very similar to the stretched 
exponential form obtained by neglecting the anisotropy 
of the segment motion in the averaging procedure. We 
propose to use the stretched exponential form also in the 
intermediate regime, where the segment motion is neither 
strictly transverse nor isotropic. This allows for a simple 



fit to the experimental data and has the advantage that 
the relevant physical mechanisms are not obscured by 
complicated numerical analysis. But, as a consequence 
of the approximations made, the numerical value of the 
prcfactor in the exponential of Eq. (2.8) should be taken 
with some precaution. Unfortunate ly, t he value obtained 
for the persistence length by Eq. ( |2.8| ) is very sensitive 
to this prefactor as well as to all the experimental pa- 
rameters entering the exponent because the persis tenc e 
length enters the expression for the decay rate, Eq. (p!s|), 

1/3 

as Lp . In addition the actual prefactor in the scaling 



law, Eq. ( |2.6[ ), for A is not known. Hence the method 
described here is not capable of producing very accurate 
absolute values for the bending modulus so far. On the 
other hand, relative changes in the stiffness, which may 
be caused by diverse chemical or physical mechanisms 
are readily detected. 



B. Initial Slope of the Dynamic Structure Factor 



Although the dynamic structure factor g{k, t) is a com- 
plicated object, and several physical assumptions and ap- 
proximations enter its explicit calculation, its initial slope 

The 



may be evaluated exactly in the case of Eq. (1.1) 
general scheme of computation may be found in Ref. |l^ ] 
(see also Ref. pl| for a summary of predictions derived 
from various semiflexible models). In the semifiexible 
case caution has to be paid to the rigid constraint of 
constant contour length. It causes an effective reduction 
of the degrees of freedom of the local quasi-free Brownian 
fiuctuations about the equilibrium configuration, which 
determine the initial slope for large scattering vectors 
and moderate chain stiffness. (For very stiff molecules 
the time regime for these quasi-free fluctuations dimin- 
ishes Q , but actin is far from this limit as can be inferred 
from a comparison of the characteristic time scales Tp, 

and 7^"'' .) For the detailed calculation of the initial decay 
rate, 



7^^ 



— log g{k,t) 



(2.10) 



for the case of a semiflexible polymer in semidilute solu- 
tion we refer the reader to Ref. Here we give the 
asymptotic result for a dilute solution and for scattering 
wave vectors k of modulus k larger than the inverse per- 
sistence length but much smaller than the inverse of 
the microscopic cutoff length a introduced in Eq. (2.5): 



'^The second condition is actually modified by entanglement. 
This may indeed be exploited in future high precision mea- 
surements, as is explained in , but is of minor importance 
here. 
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(0) 

7fe 



log ka 



(2.11) 



The deviation from ideal scaling is rather weak, hence it 
is useful to express Eq. ( 2.11 ) as a "quasi-scaling law" 

7^°-' ~ fc^^*^) with an effective dynamic exponent (Fig- 
ure 0) given by 



z{k) = 3 



6 log ka — 3 
6 log ka — 5 



(2.12) 



It is quite striking that the microscopic cutoff a appears 
here. Hydrodynamic screening does not substantially 
alter the result of Eq. ( ^.11 ) in the large wave vector 
regime, but flattens the increase of /k^ for small wave 
vectors. However, the theory is not valid for small scat- 
tering vectors, because scattering vectors smaller than 
the inverse mesh size do not resolve single filaments 
but average over several molecules. Moreover, as a con- 
cequence of local fluctuations in the mesh size there is 
scattering in the experimental data in the crossover re- 
gion k w 
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FIG. 4. Effective exponent z{k) of the "quasi-scaling" law 
for the initial decay rate of the dynamic structure factor of 
actin. 2 was predicted to be a universal number, 2 = 3, for 
all flexible polymers by classical scale invariant models. This 
was never observed experimentally. We argue that the (in- 
trinsic) semiflexibility of all real polymers is responsible for 
the discrepancy. The dashed line was computed for a semidi- 
lute solution of actin filaments [ca = 0.16 The solid 
line corresponds to the dilute case, Eq. (pd3). In the typical 
k intervall probed by QELS (5 ~ 30 Atm^ z{k) f» 2.7 for 
actin. 

Figure ^ shows a comparison of the computed initial 
decay rate 7^°'' for a solution of 0.16^ actin (dashed 
line) with previously measu red data |£^,|l9|. A simple 
least square fit of Eq. ( 2.11 ) to the data gives a = 5.4 
nm. a represents an effective hydrodynamic diameter of 
the filament to be compared with two times the cross sec- 
tional radius of gyration r^^ . The remarkable agreement 
of a with the value of 2r^ = 5.16 ± 0.3 nm determined 



by different methods [ p4[|25| ] strongly supports the above 
ideas. 
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FIG. 5. Correction to the classical prediction 7^^^' ~ 
for the initial decay rate of the dynamic structure factor. 
The theoretical predictions for dilute solutions (solid line) and 
semidilute solutions (dashed line) are compared with exper- 
imental data of Schmidt |^,^. Also included is the pre- 
diction for gaussian chains from Ref. [^. The experimental 
data and the dashed line both correspond to the same actin 
concentration Ca = 0.16 The positions of the theoretical 
curves depend on th e effe ctive hydrodynamic diameter a of 
the filament [see Eq. 2.11|. It was used as fit parameter and 
found to be a = 5.4 nm. 




FIG. 6. A real polymer is not a fractal (i.e. not a gaussian 
chain, as often assumed for computational convenience) but 
rodlike on small scales. This affects the dynamical properties 
of polymers. 

Finally we would like to note that in classical neutron 
scattering experiments with more or less flexible syn- 
thetic polymers the ratio of wavelength to persistence 
length \/Lp is not very different from what is encoun- 
tered in light scattering from large biomolecules. So one 
expects similar results in both cases. It is well known 
p3| , p^ that there is sometimes poor agreement between 
the experiments with synthetic polymers and the classi- 
cal theoretical predictions for the initial slope of g{k,t) 
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derived from scale invariant models. The excellent agree- 
ment achieved now with actin networks provides strong 
evidence that the usual scale invariant polymer models 
are not capable of describing quantitatively the dynamic 
properties of real polymers, which, after all, are semiflex- 
ible at heart (Figure ||) . 



IV. RESULTS 

We now turn to a discussion of the experimental data 
and their analysis in terms of the theory described above. 



Stretched Exponential 



III. MATERIALS AND METHODS 



Actin was prepared from rabbit muscles accord- 
ing to Pardee and Spudich with an additional 
gel filtration step as suggested by MacLean-Fletcher 
and Pollard using a Sephacryl S-300 HR column. 
Tropomyosin/Troponin were prepared from the residue 
of rabbit muscle acetone powder left after the actin ex- 
traction [^ , and separated into tropomyosin and tro- 
ponin by hydroxyl apatite column chromatography [ p9[ . 
The purity of the proteins was checked by SDS polyacry- 
lamide gel electrophoresis ||3^ stained with Commassie 
Blue, and estimated to be of at least 95% purity. Actin 
was tested for its ability to polymerize by low shear vis- 
cometry with the falling ball capillary apparatus as de- 
scribed by Pollard and Cooper and by fluorescence 
increase of 5% NBD-labeled actin Functionality 
of tropomyosin/troponin in the presence and absence of 
Ca++ was characterized by an actin binding test: they 
were added to F-actin in varying concentration ratios, 
centrifuged at 100 000 g for 1 h, and analyzed by SDS 
gel electrophoresis. 

Actin was stored in a buffer containing 2 mM im- 
idazole, 0.2 mM ATP, 0.2 mM DTT, 0.2 mM CaCla 
and 0.05 vol.% NaNa. For the polymerization of actin, 
a buffer with 2 mM imidazole, 0.5 mM ATP, 2 mM 
MgCla, 100 niM KCl, 0.2 mM DTT and 0.2 mM CaClz 
was used. The buffers where adjusted to a pH of 7.4. 
For the experiment without Ca++ the CaCl2 was left 
out and 1 mM EGTA was added. The molar ratio of 
actin:tropomyosin:troponin was 7:1:1. 

For QELS measurements, all solutions were freed from 
dust by sterile filtration, mixed, and stored overnight 
at 4° C to achieve an equilibrium polymerisation state. 
Samples to be inspected by electron microscopy were ad- 
sorbed on glow discharged carbon coated formvar films 
on copper grids, and negatively stained with uranyl ac- 
etate. 

The experimental setup for QELS has been described 
in detail previously [|l9|,^ . We used the correlator ALV 
3000 (ALV Langen) with 1024 linear channels to calcu- 
late the dynamic structure factor. The light source was 
an Innova 70-4 argon-ion laser from Coherent with 200 
mW for the 488 nm line. 



Figure show s a fit of the theoretical dynamic struc- 
ture factor, Eq. (2.S), to experimental data for a scatter- 
ing angle of 90° corresponding to A; = 24.2/im^^. Note 
that there is only one free parameter, 7^. An excellent fit 
of the experimental curves is obtained in the time domain 
10~^ ~ 10~^ s, for which the condition, Eq. (^), is ap- 
proximately fulfilled (cf. Footnote 1), whereas a simple 
exponential decay is clearly ruled out. Hence the theory 
is well suited to interpret our data within the present 
experimental accuracy. 




typical exp. result (actin 0,4 mg/ml) 
- best tit with stretclied exponentiai 
best tit with expohehtia] 



decay time in s 

FIG. 7. Fit of theoretical dynamic structure factor Eq. 
(2.5) to experimental data for k = 24.2 fim~^ . Clearly Eq. 
(2.5) describes very well the experimental situation in the 
time interval 10~^ ~ 10~^ s, wheras the simple exponential 
fit is ruled out. 



Bending Modulus 



As pointed out in Section [I A the decay rate 7fc, Eq. 
T9I), is determined by the bending modulus k and also 
by the screening length A. In order to check the effect 
of the screening length on the decay rate, we measured 
the dynamic structure factor for various actin concen- 
trations Ca- The screening length was taken to be equal 
to the mesh size ^, which was found previously p9| ] to 
obey the scaling law_| [/im] = 0.35-^Ca [^]- Th e re sults 



[IB. the 



are shown in Figure ^ As discussed in Section 
theory only applies to scattering vectors large enough 
to resolve single filaments. Scattering vectors k smaller 
than the inverse mesh size average over several fil- 
aments. Hence, the decay of the structure factor obeys 
Eq. (2.9) only for large scattering vectors k. The onset 
of the deviation may be taken as a lower bound for the 
mesh size. 

Although the concentration was increased by a factor 
of 8 the derived values for the bending modulus k are 
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fairly consistent. They agree within a standard devia- 
tion of 35% for the entire range of scattering vectors and 
concentrations. However, for a single scattering vector 
and fixed concentration the deviations are considerably 
smaller, e.g. for k = 24.2 fim-^ (90°) and 0.4 ^ the data 
are reproducible withi n 5% . Considering the discussion 
at the end of Section II A and the dependence of t he 
value of the bending modulus k derived from Eq. |2.§| on 
the hydrodynamic aspect ratio A/a - and hence on our 
choice of the screening length A - we are presently not 
able to determine an accurate absolute value for k. (An 
average value of 9.5 • 10~^^ Jm was obtained for k.) On 
the other hand, relative changes in the stiffness can be 
resolved rather precisely. 

Methods to determine A experimentally are discussed 
m Section 0. 




- actin 0.1 mg/ml 

- actin 0.2 mg/ml 

- actin 0.4 mg/mi 

- actin 0.8 mg/mi 



— r- 

25 



— r- 

30 



wave vector k [|am ] 



FIG. 8. Values for the bending modulus k obtained for var- 
ious actin concentrations Ca in the semidilute regime using Eq. 
( ^.^ ) and A = 5. Only large wave vectors k 2> resolve 
single filaments and are accessible to our theory. 



measurements at 10°C are summarized in Table |[ 



scattering vector k [fim ^] 


17.1 


24.2 


29.6 


K for actin [10~^' Jm] 


9.9 


8.7 


8.6 


K for actin + Tm/Tn -Ca++ 


11.4 


12.9 


14.7 


K for actin -I- Tm/Tn -|-Ca++ 


12.1 


15.9 


13.9 



TABLE I. Summary of values obtained for the bending 
modulus K [Jm] of actin in the presence of Tm/Tn with and 
without Ca^"'" at different scattering angles in comparison to 
pure actin. 



To exclude preparation artifacts and in order to check 
whether Tm/Tn has some effect on the mesh size, which 
would in turn affect the hydrodynamic screening length 
A and thus the derived value of k, the actin network was 
examined by negative-staining EM for all samples. We 
could not observe an effect of Tm/Tn on the network 
structure. In summary, we find that in the presence of 
Ca"''"'' Tm/Tn causes an increase of the bending modulus 
of F-actin by about 50%. 
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actin+tropomyosin+troponin - Ca^^ 
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Effect of Tm/Tn Binding 



decay time in s 

FIG. 9. Effect of Tm/Tn binding with Ca++ (0.2 mM up- 
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Figure P shows the effect of the tropomyosin/troponin 
complex (Tm/Tn) on the decay of the dynamic struc- 
ture factor for actin solutions with a concentration of 
c„ = 0.3 ^ (~ 7.1 fiM) at 10° C. Since Ca++ is known to 
regulate the coupling of the Tm/Tn-complex to actin, ex- 
periments were performed with and without Ca+''". Fig- 
ure § clearly shows the appreciable decrease of the decay 
rate in the presence of Tm/Tn. The reduction of Ca++ 
appears to decrease the stiffness slightly but the effect is 
too weak to be considered significant with the present ac- 
curacy of measurement. The experiments were repeated 
several times with different actin preparations and at two 
different temperatures (10° C and 25° C) and the same 
absolute value of k as well as the same degree of stiffen- 
ing by Tm/Tn was always observed. The results of the 



V. GENERAL DISCUSSION 

The present analysis of the dynamic structure factor 
of semidilute entangled actin solutions shows that QELS 
is a reliable tool for the study of the internal dynamics 
of se mifl exible polymers. However, as may be seen from 
Eq. (2^), the interpretation of experimental data is not 
entirely straightf orwa rd. The expression for the dynamic 
structure factor (2^) is ambiguous with respect to the 



microscopic source of an observed change in the decay 
rate. Changes may be caused by variations in the intrin- 
sic stiffness of the filament as well as by variations of the 
mesh size of the network, which in turn affect the screen- 
ing length A. As experiments with a— actinin show [ p4[ , 
even a local change in the mesh size, as caused by such 
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crosslinking proteins, results in a corresponding change of 
the decay rate over the whole range of scattering vect ors. 
This is exactly what one would expect from Eq. (2.£), if 
the main contribution to the scattered light is attributed 
to the crosslinkcd clusters (seen in EMs), which contain 
most of the filaments. 

The main uncertainty in the quantitative predictions of 
the QELS method described above presently arises from 
the uncertainty in the absolute value of the hydrody- 
namic screening length A. Setting it equal to the mesh 
size seems to be reasonable, but perhaps one could do 
better. There are various methods to determine A. The 
QELS method itself as described above could be used, if k 
was known. Direct information on A should be provided 
by the autocorrelation of a labeled filament in solution, 
which may be investigated by fluorescence microscopy in 
the case of actin. Another approach exploits reptation. 
In the same manner as the transversal friction, Eq. (2.5), 
one may introduce the kinetic coefhcient of the longitu- 
dinal center of mass diffusion of semiflexible filaments in 
entagled networks C|| = 2'Kr]L/\og{K/a). 

There are also effects from the center of mass and ro- 
tational motion, which we have neglected so far. Because 
of the polydispersity of the actin solution, some filaments 
shorter than the mesh size [L < ^) that are free to ro- 
tate and diffuse are always present. Their hydrodynamic 
correlation is not disturbed, nor are they hindered by en- 
tanglement. The center of mass and rotational dynamics 
of those short filaments may to a good approximation be 
described by the theory for rigid rod molecules ||l^ (for a 
more complete treatment including bending motions see 
1^). For the time being, the accuracy of the experimen- 
tal data does not allow for a quantitative analysis of this 
small effect. 

Finally we want to compare our findings with results 
obtained by the fluorescence technique. The two methods 
can be considered complementary: They probe adjacent 
wavevector and time regimes. The fluorescence technique 
is more direct. QELS, on the other hand, avoids prob- 
lems associated with potential perturbations from fluo- 
rescence markers (for the case of actin and phalloidin see 
Ref. [Q) and gives much better statistics. For the time 
being, a significant discrepancy remains concerning the 
effect of Tm/Tn. Tm/Tn was seen to cause a softening 
of the filaments with the fluorescence technique Q . This 
contradicts the QELS data given above as well as some 
recent direct measurements |3|,^^. We do not have a 
fully convincing explanation for this discrepancy. A stiff- 
ening of actin by Tm/Tn binding seems very probable, 
since other proteins binding along actin filaments are also 
known to enhance the stiffness of the filaments, e.g. talin 
P?! , a protein involved in the membrane binding of actin 
in cells. 

However, one may also think of a more intriguing sce- 
nario in terms of torsional modes and other nonlinear 
or higher order derivative contributions to the configu- 
rational energy of the filament. No thorough theoretical 
treatment of all those possible contributions has been 



given so far, and it is not known, which of them are ac- 
tually relevant. A first principles analysis of the various 
linear modes of actin (in vacuo) has only been published 
very recently js^. Especially for F-actin local torsional 
modes and additional "groove swinging" and axial slip- 
ping motions are thought to be important because of its 
double-stranded structure. Very strong evidence for local 
torsional excitations was provided by high resolution elec- 
tron microscopy studies showing that the twist angle 
is not fixed but may fluctuate by ±10°. Torsional fluctu- 
ations of some 100 nm length have been identified. The 
presence of torsional modes may affect the dynamic struc- 
ture factor in two ways: First, they cause "crankshaft" 
motions of bent chains and thus provide an additional 
mechanism contributing to the time decay of the correla- 
tion function. Second, they could lead to a scale depen- 
dent renormalized bending stiffness of actin as was re- 
cently observed for the railway track model . Though 
we did not see any systematic ^dependence of the bend- 
ing modulus with the QELS method, such effects may 
well occur in other /cregimes. The more "exotic" mo- 
tions of actin are of great biological interest, since they 
could play an active role in the conformational changes 
of the actin-myosin-complex during ATP-cleavage. It is 
hoped that future QELS studies with enhanced precision 
and larger ranges of scattering vectors and measurement 
times will help to clarify this important point. 
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